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I. Introduction 

A. The Problem and Its Importance 

As knowledge of man's physical environment is improved, the relative lack 
of information about the seas around us becomes ever more important. We 
actually know far more about the moon's surface than we do about the waters 
beneath the surface of the ocean. 

The knowledge, however, of the condition of the water around us is con- 
tinually shown to be vital. Recreation, fisheries, transport, fresh water 
utilization; all these are threatened by various forms of pollution of the 
waters, as depicted in Fig. I.l. To cope with this situation, accurate 
measurements of various water parameters must be provided. Although gross 
indications of major ocean currents and water masses are known, there has 
been until very recently virtually no continuous source of information about 
even such elementary quantities as sound velocity, temperature, salinity and 
turbidity. Without a knowledge of these basic building blocks of any theory, 
we shall find it difficult to proceed in our conquest of the water surface of the 
earth, which comprises more than three fourths of its surface. Very recently, 
it has become possible through the methods of remote sensing to obtain some 
of the needed data. By measuring the infrared emission from the water, the 
surface temperature may be continuously determined. The water surface radiates 
i nf **ared energy very nearly like a black body (a perfect absorber), ihe inten- 
sity of this radiation at a wavelength of 10 micrometers can be related to the 
sea surface temperature. However, infrared emission of this wavelength can only 
be used for temperature measurements of the upper few hundredths of a millimeter 
of the water, since water is extremely opaque to radiation of this frequency. 

Often, the temperature of this upper skin is not really representative of 
the body temperature of a body of water. This is because the processes by whicn 
water loses its heat occur in the upper tenth of a millimeter or less immediately 
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below the surface. These processes include radiation, evaporation and con- 
duction. The result is that, unless there is violent mixing induced by a 
strong wind, the upper skin of the sea may be a degree centigrade or even 
more cooler than the water just below. ^ It is this lower temperature that 
is measured by ordinary infrareo remote sensing, so that a true picture is 
not obtained. 

In addition many problems in oceanography require data not only from 
the surface but from deeper layers in the ocean. Typical of these prob- 
lems is the question of survival of life in the neighborhood of a fossil 
or nuclear power plant cooling canal outfall. Here, for example, if the 
deeper waters remain tolerably cool, the damage is minimized, while if the 
water beco.mes too hot all the way to the bottom the result is catastrophic 


Another important piece of data is the existence of a region of rapid change 
in temperature beneath the surface. Such layers are common in the sea and 
cause problems in acoustic transmission, since they bend or reflect sound 
waves; their detection is often important, especially in military applications 
As present, the only way to measure the subsurface temperature at depth 
is to introduce thermometers directly into the water to the desired levels. 
Such methods are expensive, time-consuming, or do not provide the continuous 
picture that remote sensing could provide of the surface. It would be very 
valuable, therefore, to be able to provide remote sensing of the ocean 
parameters at various depths. The methods described here will provide con- 
tinuous information on temperature, salinity, sound velocity and backscatter 
turbidity without directly introducing devices into the ocean. Ultimately 
remote sensing may be carried out from an aircraft. Initially, however, the 
measurements will be made from a platform such as a dock or a ship on the 
water surface. Other important oceanographic parameters to be measured 
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include salinity, turbidity and the velocity of sound. These parameters are 
depicted in Fig. 1.2, with representative curves for the principal subjects of 
this report, temperature, turbidity, and salinity. Mote that each parameter is 
determined as a function of the depth, h. In addition, information on total 
depth and the sea surface may be obtained. 

These quantities are measured at present by relatively slow and tedious 
methods. The classical procedure involves the use of Mansen bottles, as shown 
schematically in Fig. 1.3. A cable is lowered with a high capacity winch, 
and the Nansen bottles are attached by a technician standing on a grating sus- 
pended over the sea. Each bottle is equipped with a thermometer, a pressure- 
sensing device, and stoppers. When all the bottles are in position, a messenger 
weight is dropped down the wire. The weight strikes a trigger on the first (top) 
bottle which closes the bottle, activates the pressure sensor and the thermometer, 
and releases a second messenger weight which falls to the next bottle, and the 
process repeats until all the bottles are closed. The string of bottles is then 
raised by the winch and each bottle removed as it passes the grating. The tem- 
perature and pressure sensors are read and the bottle is emptied into a flask 
which is stored until analysis in the laboratory (for turbidity and salinity, 
for example) can be performed. This is a time-consuming process. For the 
measurement depicted (samples each 5 meters for a depth of 40 meters), approx- 
imately one hour would be required for the station. The ship must, of course, 
be held as motionless as possible during the operation. 

A more modern and rapid method of obtaining the oceanographic data is by 
the use of temperature, pressure, salinity and turbidity sensors which measure 
these qualities electrically, transmitting the information to the surface 
continually via a suspension cable, which carries insulated electric wires in 
its core, see Fig. 1 . 4 . This method is far more raoid than the Nansen bottles 
for the same deoth as before, perhaps only 15 minutes is necessary for each 
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station, not counting the time to stop the ship. The measurements, however, 
are not simultaneous, since the sensors must be mechanically raised or lowered 
between each depth. Again, however, the ship must be held strictly stationary 
during each series of measurements. The process of slowing and stopping a 
ship is very time-consuming, in many cases taking almost as much time as the 
actual making of observations, 

B. Advantages of Our Proposed System . 

Our method, first described in detail at the Symposium: "Remote Sensing 

Applied to Energy-Related Problems", Miami, Florida, December 1974 and 
reproduced in the Proceedings, is far more rapid. It utilizes a laser beam 
to probe the sea, and the measurement is essentially carried out at the 
velocity of light in water, 0.75 c. As is shown in Figure 1.5, in this 
method a laser beam is projected downward into the water. Return light is 
received, and the information contained is computer-processed to obtain the 
speed of sound, temperature, turbidity, and salinity, all as a function of 
dapth. The ship need not stop at all, even during measurements. This makes 
the new method very much faster than any previously used procedure. In five 
minutes, for example, a ship steaming at a steady speed of 7 knots can make 
■|50 measurements (one each 2 seconds), with a separation of about 7 meters 
between each observation. Such a series with the sensors would require at least 
15 minutes per sounding, assuming 10 meter depth, or about 35 hours, and with 
bottles considerably more! 

In addition, the proposed method, as we will show, can be adapted to making 
airborne measurements. In this case the saving in time is further extended 
so that rapid simultaneous surveys will be possible, which at present are far 
beyond the range of practicability. 

As an example of the kind of problem that could be studied with 
a system of this kind would be the three-dimensional behavior of an estuary 
under conditions of tidal flow. Tides change each six hours in typical 'J.S. 
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East Coast estuaries, and maximum currents and slacK conditions often last for 
a very short time^ (a few minutes in many cases). A three-dimensional grid 
of aircraft observations of temperature salinity and turbidity can then be 
compared to three-dimensional theoretical models such as that developed in 
NASA Project NASlO-8740, "The Application of Remote Sensing to Detecting 
Thermal Pollution", University of Miami, Coral Gables, Florida, 1974 to date. 

In this way our knowledge of oceanographic parameters in a given dynamic 
situdtion will b6 incr^stssd msiny-folcl. It would, thsrstors, sesm thdt ^ Qusn- 
turn jump in understanding of the dynamics of :he sea could result from the 
application of this method. 

C. Summary of Feasibility Study . 

Any proposed method, no matter how attractive in theory, must be studied 
in detail before too much time and money is lavished upon it in order to see 
if it is feasible at all. Accordingly a Feasibility Study was made by our 
Laboratory, which was reported in the NASA CR-139184 document^ dated January 
1975 mentioned above. 

This study examined whether there was enough intensity in the returned 
scattering from the sea to make measurements of the parameters, turbidity, 
salinity, and temperature. Using information theory, this data was then 
expressed in terms of accuracy obtainable in a given time of observation. 

To obtain information from layers beneath the surface, it is necessary 
to use radiation for which the sea is transparent. The band of electromagnetic 
radiation which penetrates jest into seawater lies between 0.48 and 0.55 
micrometers, or between 4800 and 5500 Angstroms*. It is, therefore, proposed 


*0ne Angstrom equals 10"^'^ meter. 
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to project a beam of light this wavelength into the sea and to use it 
to probe the deeper layers for information on oceanographic parameters, par- 
ticularly temperature, salinity and turbidity. Either of two methods can 
be used. A light source constant in time can be projected downward while 
a collimated detector system some distance away is used to vary the position 
in depth of the scattering received. This is the bistatic system which, 
being the simplest, will be the first investigated experimentally. An 
alternate system employs col inear source and receiver beams; here the depth 
information is obtained by means of the time delay between projected and 
received signals, analogous to radar and lidar. For this purpose, rapid 
modulation of a continuous light or short pulses of light (of the order of a 
few nanoseconds) are necessary in order to provide meaningful depth infor- 
mation. A moderately powerful light source is also required so that measure- 
ments may be made to as great a depth as possible. Finally, the temperature 
and salinity information is contained in small changes in wavelength of the 
light returned from the sea, so that the light source used must have a narrow 
wavelength spread. All these requirements can be satisfied by the use of a 
laser. The laser will be mounted on a ship or aircraft and will project its 
light downward into the sea. Returning light scattered in a backward direc- 
tion by the sea will be detected and analyzed with respect to time, intensity, 
polarization and wavelength. This information will in turn be recorded and 
processed to yield the desired parameters as a function of depth. The system 
is shown schematically in Fig. 1.6. 

In order to measure the parameters the evaluation of light scattered by 

a number of processes is necessary. 

1 ^ Scattered light from a laser beam : 

a) Tyndall and Mia scattering: Here the wavelength of the incident 

light is unchanged. The intensity is, in general, proportional to the density o 
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particles suspended in the water. These forms of scattering give rise to 
cloudiness or turbidity which is the phenomenon that interferes with the 
clarity of the water. The amount of scattering depends on angle, but for 
most situations of interest a single measurement at a particular angle can 
be used. Here we propose a normalization function to correlate measurements 

made at different angles. 

b) Rayleigh scattering: This occurs in clear water and is analogous 

to the molecular scattering in the atmosphere that gives rise to the blue 
color of the sky. The intensity of Rayleigh scattering is usually much less 
than Tyndall scattering, and its intensity is a known function of tempera- 
ture and can be predicted. In general. Rayleigh scattering will be of no 
direct use to us. but because it can be predicted and is very faint, it will 

not interfere with our measurements. 

c) Brillouiti scattering: Photons in the Incident light react with 

phonons in the liquid to produce wavelength-shifted light, where the shift 
is proportional to the velocity of sound in the liquid. The shifted wave- 

length'* is given by: 


X = Xq (1 ± 2n ^ sin ) 


where is the incident wavelength, n is the index of refraction of the 
water, v^ is the sound velocity, c is the velocity of light, and e is the 
scattering angle. This means that there are two shifted spectral lines with 
scattering proportional to v^. the sound velocity. For backscattering. the 
angle between incident and scattered radiation is .. The resulting wavelength 
shifts are about 0.08 A, which can oe measured using a Fabry-Perot inter- 
ferometer. This measurement yields the sound velocity in the sea, v^. 

The Fabry-Perot is a spectroscopic instrument with extren-.ely high re- 
solving power and sensitivity. In addition it can be prepared^ to have several 
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channels operating siinultaneously. Thus the measurement of the Brillouin 
peaks can be carried out rapidly enough to preserve the depth information 
provided by the time delay, which amounts to 6 nanoseconds per meter. 

d) Raman scattering: The Raman scattering differs from the com- 

ponents considered so far in that the wavelength of the scattered light dif- 
fers much more from that of the incident. For this reason, the spectroscopic 
problem is much simpler, and interference filters can be used for its measure- 
ment. Salinity information can be obtained either from depolarization of 
water Raman lines, or the intensity of certain ion lines, principally the 

sulphate ion. 

In the Feasibility Report it was shown that a moderately powerful laser 
operating in the blue region of the spectrum (440 to 490 nanometers wavelength) 
can, for a relatively clear sea, be used to provide data on the parameters men- 
tioned above as a function of depth. Brillouin, Raman, and Tyndall scattering 
from the laser beam are all to be utilized, as summarized below. The process 
for each parameter is described in turn, 

2. Sound Velocity 

The sound velocity, v^, in the sea is well known, having been of military 
as well as scientific interest for many years. It depends on salinity, depth 
and especially temperature. It turns out that the sound velocity in 
liquids can be measured using a beam of light as a probe, and that light from 
a laser is particularly well suited to the purpose. If such light is passed 
through water, it is found that two additional wavelengths appear in the scat- 
tered light. These additional wavelengths carry the information on the speed 
of sound; their spacing is proportional to v^, the sound velocity. For back- 
scattering, the wavelength shifts are about 0.08 A, which can be measured using 
a Faory-Perot or a Mach-Zehnder interferometer. This measurement yields the 
sound velocity in the sea, v^. 
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3. Salinity and Temperature 

Measurements of the salinity are to be carried out by using the Raman 
effect. In the Raman effect, unlike the Srillouin effect, the wavelength of 


Che backscattered light is considerably altered by the interaction of the 
light with energy levels of the salt and water molecules. Since the population 
of these levels is temperature dependent, the characteristics of the Raman 
spectrum depend on the water temperature and the salinity. It turns out that 
both the wavelength and the polarization of the Raman light are temperature 
and salinity dependent. A recent study® by Chang and Young has shown that if 
the laser light is circularly polarized, the Raman scattered light will show 
depolarization which is dependent on salinity and temperature. Tnis dependence 
is different from that of the Srillouin effect mentioned above. A combination 
of the two observations yields the temperature and salinity seoarately. 

A more direct method involves the 30^ ion. Due to its internal structure 
this ion gives rise to a vibration-rotation Raman effect which is favorable 
in wavelength to transmission by the sea. Except in very unusual circum- 
stances, the 30^ concentration and those of the principal constituents of 
the salinity of the sea, Na and Cl, are strictly proportional . This means 
that a measure of the 30^ ion will suffice to measure the salinity, ihere 
seems no reason why this cannot be done with sufficient accuracy to reduce our 
temperature measurements. An exper^nental study of this approach in seawater 
is planned. The salinity measurements are then combined with those of sound 
velocity to obtain the temperature. 

4. Turbi di ty 


In the literature there are at least tan definitions 


of turoidity. 


We 
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choose it to be the lack of clarity of water resulting from suspended 
particles which causes light to be scattered at random. A laser beam is 
particularly well suited to make this kind of measurement, since the quantity 
of backscattered light from a laser beam is a direct measure of turbidity. 

This component of scattered light exhibits no change in wavelength from the 
incident light and is called Tyndall scattering. The ratio of the intensity 
of Tyndall scattering to incident light intensity gives the turbidity. In 
previously used methods a difficulty has always been the incident light in» 
tensity measurement. The intensity of the Brillouin scattering corresponding 
to a given incident intensity is a known function of the water parameters and 
will give a direct measure of the intensity of the laser light at a particular 
layer of the sea. This is then compared to the Tyndall component, and the 
turbidity of the various depths can be directly deduced from their ratio. 

Thus, the measurement of turbidity can now be made much more accurately than 
before, since a ratio is used where both quantities come from the same water 
sample. 

D. Feasib'jlity Criteria 

1) Seaborne Operation. 

The acid test of any measurement is the accuracy obtainable in a given time 
of observation. The mechanisms of Tyndall, Rayleigh, Brillouin and Raman 
scattering are sufficiently well known that a calculation can be performed 
on the beam of scattered light from the sea when the optical parameters 
(size of lenses, for example) are known, and the intensity of the incident laser 
light is given. 

We have assumed the use of a relatively modest laser {0.1 joule), an 
average optical system, and determined the accuracy of each measurement by 
means of a computer simulation. First the intensity of scattering was determined. 
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This calculation was checked against a 3rillouin spectrum in the Laboratory, 
as Shown in Fig. III. 2 (p. 32). Then using information theory and the 
characteristics of typical interferometers of the type to be used, the in- 
formation content of signals from the sea was evaluated. Fig. UI.4 (p. 34) 
shows a computer simulation with a separation of the 3riHouin peaks of 
0.28 cm*^ and a turbidity of 0.235. 

The main configuration considered was a three channel Fabry-Perot inter- 
ferometer. shown in Fig. VI. 6 (p.68). A laser sends its light by a right- 
angle mirror into the sea. It is backscattered and collected by the collecting 
lens. Rendered parallel by a negative lens the light passes through a Fabry- 
Perot interferometer and a focussing lens and falls on a Fafmr image divider. 
The three channels are separated here and are directed to separate onoto- 
multi pliers as shown. 


a) 


The parameters chosen were: 
Power 

Wavelength 
Line Width 
Pulse Length 


LASER 


100 mj/pulse 

O 9 

variable 5,300 A - 4,000 A 
0.05 cm"'* 

4 nsec 


b) 


GEOMETRY 

Height Above Surface 
Sea Surface Factor 
Area of Collecting Lens 
Efficiency of Instrument 0.1 


0 (Shipborne) 

1 

300 cm^ 


A typical result is shown in Fig. VII. 1 (p.32). Here tne accuracy of temcerature 
measurements is shown at various ceoths in clear seawater as a 'unction or wave- 
length. It is seen that the most favorable wavelength lies in tne neignoorncod 
of 450 nanometers ana at 50 meters an accuracy of about 0.2 degrees ■. s 
forecast. 
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One of the most important questions pertaining to an optical system of 
this sort is the effect of turbidity in the water on the depth which can be 
reached. This can be answered by the computer model and a typical result is 
shown in Fig. VII. 3 (p. 84). Here two Fabry-Perot interferometers, one with high 
finesse F and one lower are compared. Accuracy in degrees centigrade obtainable 
at 20 meters is plotted against turbidity. We note that the effect of turbidity 
is relatively small, indicating that clarity of the water is not as vital as we 
had thought, an encouraging result. 

2) Airborne Operation 

For the airborne experiment, the height above the surface Csee Parameters, 
part b, above) may be given a value, for example, of 200 meters. Here the 
effects of sea-surface also enter. However, the airplane bathymetry experi- 
ments of the Wallops Island group in the Key West area have indicated that, 
except for very strong wind conditions, the sea surface does not appreciably 
interfere’ with the bathymetry accuracy. Our own calculations indicate the 
same th-<ng. The feasibility of airplane measurements is certainly indicated. 

3) Bi static Experiment 

The airborne measurements require a coaxial laser and scattering team, 
so that the effect of the non-flat sea surrace is eliminated. This, in turn, 
means that depth information is to be obtained from time-of-flight data. Such 
measurements require fast pulses from the laser as well as very short gated 
receiving systems. These lasers are expensive to buy and probably cannot be 
obtained surplus. 

It happens, however, that for a seaborne experiment another far less 
expensive method can be used for depth information. This is to use two 
positions, one for sanding and one for receiving, and to vary the angle be- 
tween them for depth measurement. This is the principle of the ordinarv 
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rangefinder. Such a configuration is shown in Fig. VIII. 1 Cp- 881- The 
depth of the shaded area from which the measurement comes can be deduced 
from the angles of the two beams by trigonometry. Our first experiments in 
actual outdoor salt water will be made using this configuration. This is 
because due to a restricted budget, we are relying mostly on equipment either 
on hand in the Laboratory or borrowed Oike the laser). Since this method 
measures the Brillouin, Raman and Tyndall Scattering as a function of depth 
just as the final system will do, it provides an extremely valuable check 
on our calculations, and so is a logical first step in our Experimental 
Program. When we have actual Brillouin spectrum results from a natural 
body of water outside the Laboratory, our level of confidence will be greatly 
increased. 
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II. Oelineatiofi of the Quantities 
A. Temperature 

The object of our measurement here is the thermodynamical bulk 
temperature of the water, the quantity normally measured when a thermometer 
is introduced into the water. This temperature is required in considera- 
tions of water circulation and biota distribution and can be used as a 
labelling parameter of water masses in certain circulation studies. 

In our measurements, -.he temperature is inferred from the speed of 
sound, and, assuming sufficient knowledge of the other parameters pressure 
and salinity, this inference is very precise. An expression we have used 
for the speed of sound is exact Jo better than 0.5 m/s®, and with the linear 
temperature variation of about 3 m/s per ®C, the precision is already better 
than 0.2 ”0. Using suitable calibration, we do not expect any observable 
error in this inference; the temperature will be computed with the precision 
with wh'^ch the speed of sound, salinity and pressure are known. 

The Brillouin spectrum is a very direct method of measuring the speed 
of sound: the wavelength of the laser pulse is shifted in the scattering 

volume by an amount directly proportional to the speed of sound at this 
location; this modified light is then observed without any possibility of 
further change, since a second scattering would reduce the intensity to 
nothing. For comparison, the only presently employed method of temperature 
remote sensing by thermal infrared intensity measurement is sensitive only 
to the observed skin temperature. There are indications^’® that this skin 
temperature may be lower than the bulk temperature by l^C and even more, 
depending on relative humidity and wind conditions. In addition, the in- 
frared measurements have to be corrected for the absorption in the atmosphere, 
which varies with path length and humidity. 
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The measurement of temperature using the Brill ouin spectrum^on the 
other hand, is only limited by the power of the laser and the ability to 
measure salinity. 

B. Turbidity 

Until very recently, the quality of our waters in lakes, streams and 
near the coast at sea has been deteriorating steadily. Dredging operations, 
discharge of industrial and municipal wastes, and agricultural usage all have 
contributed to a decrease in quality of those very neighboring waters most 
useful to us for recreation and for fisheries. Since the general realization 
a few years ago that we must preserve our water heritage and limit our 
despoliation of this important natural resource, there seems to have been 
a dramatic turn of events. In some cases tie rate of deterioration has 
apparently been slowed, and even reversed. Controls on activities which 
may adversely affect water quality have multiplied, and considerable expen- 
ditures have been made or projected to improve water quality. Fisherman 
have reported improvements in the Great Lakes, and coastal waters such as 
Biscayne Bay seem somewhat clearer than a few years ago. But how is water 
quality to be measured? How are we to know whether we are winning the battle 

whether the expenditures are fully effective. 

Of the measurements of water quality, that of turbidity or lack of 
clarity is probably the most useful. This is not only because measurements 
can be made optically, and therefore quickly without complicated chemical 
analysis, but also because it is exactly the lack of transparency which 
turbidity implies that is perhaps the most direct and obvious measure of 
water deterioration. 

The literature discloses^® at least nine distinct definitions of 
turbidity or lack of clarity. Here we shall consider that turbidity causes 
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the cloudiness in waterwhich obscures objects on the bottom. We shall not 
be concerned with absorption, which limits the transparency of water but 
does not reduce the definition or contrast of what we see. The cloudiness 
that we wish to measure is caused by suspended small particles in the water. 
When light encounters these particles, it is scattered in all directions. 

The result is a decrease in the clarity and contrast of an image. 

Consider a beam of monochromatic light incident on a sample of water. 

It is assumed to have wavelength X and irradiance In an infinitesimal 

volume dxdydz, a certain amount of light, dl, is scattered in a solid 
angle dn. This amount is proportional to the incident intensity U, to 
the solid angle dfi, and to the volume dxdydz, and is a function of angle 
and wavelength: 

dl = a(9,x) lo dxdydz 

Here a(9,x) is the scattering coefficient which is to be measured. In 
those cases where the particles doing the scattering are small, of diameter 
d, the same order or smaller that the wavelength x of the light, the angular 
dependence is very complicated. In most of the turbidity which affects 
water quality, however, the particles are considerably larger than wave- 
length, d » X, and the variations with angle are simple while the wavelength 

dependence is negligible. 

Measurements of turbidity are thus equivalent to the measurement of 
scattering, which is dependent on angle. A type of scattering meter, the 
nephelometer , is adapted to measurements at 90”. Other methods include 
the taking of water samples, the Secchi disk, and measurement of the 
intensity of a light beam after passage through a definite length of water 
compared to its original intensity. They may be classified as follows: 
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1) transmissomster 

measure either diffuse or true extinction coefficients 


either in situ or in vitro. 

2) scatterometer (nepheiometer) 

measure the scattering properties of the sample, usually 
but not always at 90° in vitro. 

3) combination of the preceding two methods 

a) contrast meter; - Jackson candle turbidimeter 

- Secchi disk 

- rotating black and white disk 

b) ratio of transmitted beam and scattered beam. 

In order to be able to make even the most approximate comparison of 
measurements by these methods, some kind of standard angle dependence must 
be assumed. We found that instead of making each of the comparisons with 
ad hoc assumptions, it is much simpler and more consistent to refer all 
these measurements to a new unit based on the definition of turbidity; 
"Turbidity corresponds to the amount of suspf ded 
matter in a sample of water, as ascertained by op- 
tical observations. The units are required to be 
proportional to dilution, viz; dilution by half 
of a sample will reduce its turbidity by half. 

We define then a sample of turbidity 1 as having a scattering coefficient 
<j(e,x) equal to the turbidity normalization function 3(5) shown in Fig. II. i 
This function was obtained as the simplest closed analytical form that 
approximates the scattering angular dependence of most bodies of water. 

Its intensity corresponds roughly to the transition between clean and 
dirty water. Some technical data of this normalization function are 
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are given below, and all the turbidities discussed in this and the 
following reports will be given in this unit. 

The normalization turbidity function 3(9) : 

1) defines a sample of turbidity 1 and has an angular dependence 
chosen to approximate the scattering of most bodies of water 
(for instance, 

, 0 ^ . 3 

4it a 

while this quantity varies from 2.5 to 3.5 in oceans and lakes.) 

2) when multiplied by 16.77197 gives an angular function in cur- 
bidity normalized to 4 tt. 

3) can be understood as corresponding approximately^^ to about 

1.2 X 10® particle liter, of index 1.6, witn diameters 
above 2 urn distributed in Junge distribution slope 1 (assuming 
density of 2.5, this gives about .6mg/liter). 

In the ocean, the actual particulate matter giving this 
scattering would probably be closer^^ to 5 x 10'-® particles/ 
liter, of index (1 .4 -t- 0.01 i ) . distributed between the diameter 
of 0.08 and lOum by a Junge distribution of slope 3. This 
results in an absorption coefficient a * 5.2 x 10*^ that we 
use as the normalized absorption coefficient o^ 9 sample of 
turbidity 1. 

4) the total scattering coefficients of a sample of turbidity t 
is 1^ X 0.59623; with the normalized absorption i of the par- 
ticles, this corresponds to a half intensity propagation 
1 ength HIL » cm. 
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C. Sal inity 

The presence of salt in sea water results in appreciable changes in 
the physical properties of the water. The simplest expressions for the speed 
of sound, the index of refraction, and the density must have at least one 
salinity parameter. Fortunately, since most of the variations of salinity 
are due to rain or evaporation, the ratio of the different salt constituents 
stays about constant, and only one salinity parameter is quite sufficient 
for most uses: most physical properties of sea water can be written in a 

T p S system, i.e. as a function of only temperature, pressure and salinity. 

The salinity parameter can be expressed in many ways, but parts par 
thousand of salt has been a fairly standard one for a long time. It has the 
advantage that it can be understood as gram of salt per kilogram of sea 
water in thermodynamical computations. These are the units we use in this 
report. A normal value for salinity in the ocean is 35 gm/kgm and is 
distributed as^**: 

19.353 Chlorine 
10.76 Sodium 
1.292 Magnesium 
2.712 Sulfate 
0.^112 Calcium 
0.399 Potassium 

The rest is Bicarbonate and some trace elements. 

0. Speed of Sound 

The Brillouin spectrum corresponds to two peaks shifted in wave- 
number by an amount equal to 

2 n w , where n is the index of refraction 

Ly sin ie is the speed of sound 

c 0 2 is the laser wavenumber and 

is the scattering angle. 
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Since the wavelength of a laser is well fixed and the sin ^ ^ variation 
is negligible in backscattering (a » 180), the separation of the peaks 
gives a very precise measure of n v^, or more exactly the index of refraction 
at the laser wavelength multiplied by the speed of sound at that frequency: 

2n V V sin i 9 . This is a very well defined physical quantity, which is 
so 2 

a local bulk property of the water at the scattering point. Although this 
water characteristic can be referred to the T p S system, its measure is very 
direct and, since it is never influenced by the surroundings, it may be 
advantageous to consider it independently in certain studies. We expect 
these kind of measurements to be very useful in studies of stratification 
of water bodies, to obtain a fast three dimensional picture of the plume 
from a river or cooling water outlet, and as a labelling parameter in cir- 
culation studies. 
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III. Turbidity 

A. Theory 

A computer model of the Brillouin Spectrum of water based on the theory 
reviewed in Chapter V of this report was developed during this phase of the 
contract. A brief review of the Tyndall or particle scattering information 
used is presented in this section and is shown to be a means of measuring 
turbidity when coupled with the intensity of the Brillouin spectrum. 

Tyndall scattering was mentioned^ as typically the brightest of the de- 
tected returning radiation. Its central wavelength and spectral width are 
essentially the same as those of the laser. The angular dependence of the 
scattered intensity has been studied by others in many bodies of water. The 
combining of results of several of these experiments has led to the turbidity 
normalization function shown in Fig. II. 1 (p. 24) which matches tne data 
over a wide range of variation. 

The separation of the Brillouin peaks contains information about the 
water sample temperature and is independent of the absolute intensity of 
the laser. Since the temperature of the water is calculated, it follows that 
the intensity of the Brillouin effect carries the information needed to deter- 
mine the intensity of the probing laser light that interacts with the vo.ume 
of water being studied. Using the turbidity normalization function mentioned 
above, the normalized Tyndall peak intensity can be directly compared to the 
Brillouin effect intensity and an accurate and repeatable measure of turbidity 

obtained. 

This method of measuring turbidity has a significant advantage over 
previous methods because the standard for comparison (Brillouin) 
originates in the same volume of water as the effect under study (lyndall). 
Both the Brillouin and Tyndall scattered light reaches the receiver after 
transversing the same optical path so no assumptions need be made about the 
medium between the receiver and the volume under study. 
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B. Experimental 

During the feasibility study of this project the Brillouin spectrum 
was observed using a helium-neon laser and a Fabry-Perot interferometer. In 
this phase of the project the computer model was used to simulate the 
laboratory conditions as reported earlier. The purpose was to obtain "a 
posteriori" the turbidity of our water samples by matching the laboratory 
data with the computer simulation. 

1 . laboratory Work 

The light source used was a 3PECTRA-PHYSICS Model 124 helium-neon laser 
(6328 A) with an output power of 10 mW. The water samples were contained in a 
spherical glass flask and the purity of the sample was varied. The non-specific 
description of purity was then replaced by our definition cf turbidity so 
in essence, scans of the Brillr n spectrum were made on samples having dif- 
ferent values of turbidity. Most of the scans were performed at a scattering 
angle of 90®. A lens system relayed the scattered light to a Fabry-Perot 

interferometer as shown in Fig. II 1. 1. 

The Fabry-Perot interferometer was in a pressurized housing with a 1 cm 
spacer giving a free spectral range of 0.5 cm“^ (0.20A at 6323 A). yJith this 
spacer the Brillouin scattered light fell approximately 0.3 free spectral range 
from the Tyndall central peak. Freon 12 (€ 2012 ^ 2 ^ 

Fabry-Perot with a change in pressure of about 0.3 Ibs/in^ corresponding to 
a complete free spectral range of the interferometer. 

A lens imaged the fringe pattern on a plate with a centered small aperture 
followed by a photomultiplier. When the laser light was put through the 
optical system, the finesse of the whole system was measured and found to be 
about 5. 
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A typical output from this system including the Brillouin spectrum is 
shown in Fig. III.2. The scattering angle was 90® and the water was considered 
to be relatively clean. The limitation of available finesse Cusing off-the- 
shelf components only) did not permit better isolation of the central peak. 

It is moreover not necessary to clearly isolate the Brillouin peaks to obtain 
an accurate measure of turbidity. A signal-to-noise calculation in chapter '/II 
of this report demonstrates this effect in detail. Typical results indicate 
that the method easily gives turbidities with a precision of a few per cent, 
which is much better than any previous method. 

2. Computer Simulation 

The computer model of the Brillouin spectrum was developed using the 
theory presented in Chapter V of this report. The temperature, salinity, 
and pressure of the water are input data. From these specified quantities 
other properties of the sample are derived. These are density, index of 
refraction, speed of sound, partial derivatives of the index of refraction 
(with respect to temperature, pressure, and salinity), volume expansion coef- 
ficient, and isothermal compressibility. 

Additional information is then added, giving the angle of scatter of 
the receiver and the wavelength of the laser. The scattering coefficients, 
spectral widths, and spectral shifts for each type of scattering are then 
calculated (isobaric, salinity, Brillouin, and anisotropic). A Tyndall 
scattering coefficient can then be assumed and the complete simulated 
spectrum displayed. 

The simulated spectrum can be put through a theoretical Fabry-Perot 
interferometer. This can be done after a finesse and free spectral range 
is specified. Figs. III. 3, III.4, and III. 5. show a sequence of modelings 
as the scattering coefficient rp is varied. The pressure, temperature, and 
salinity were assumed to be 1.0 atm. 20®C,and 0.0 parts/thousand respectively . 
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FIGURE 111.4 



FIGURE III. 



The scattering angle was 90^ and the laser wavelength was 6328 A. The Fabry- 
Perot was assigned a free spectral range of 0.5 cm ^ and a finesse of 3. 

The closest match (relative height and half-width) is Fig. III.A with 
a turbidity of 0.285. 
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IV. The Problem of Determination of Sa1irv=ty_ 
A. The Need for the Measurement 


Brillouin spectroscopy can provide a measurement of the speed of sound 
in the sea. This is discussed in more detail in other sections of this report. 

The speed of sound is given as a function of temperature, salinity, and 
pressure (depth) by many authors. The depth for each value of the speed of 
sound can be found by measuring the time delay in the arrival of the back- 
scattered light from the sea following the arrival of a reflection from the 
surface; or in a bistatic system one can provide depth information by the use 
of trigonometry. 

The temperature and salinity information cannot be separated by the 
speed of sound measurement alone, and so an independent method for measuring 
one or both quantities must be achieved, except in special cases described 
below, where simplifying assumptions can be made. 

1. In fresh water the salinity may be assumed to be zero. The speed 
of sound and depth measurements are then sufficient to provide a measurement 
of water temperature. 

2. In the ocean the salinity near the surface changes gradually over 
distances of hundreds of kilometers. Figure IV. 1 depicts the mean surface 
salinity of the world's oceans. 1“ The surface layers are often well mixed 
and therefore the salinity may be assumed to be constant. A mean value for 
salinity may be chosen with which to calculate temperature from speed of sound 
measurements, allowing for some uncertainty of results. 

3. Some water bodies of interest may not fit these special cases. 

For example, brackish water in estuaries generally exhibits large vertical 
and horizontal gradients of both temperature and salinity, as does the sea 
in regions of high precipitation or evaporation. In these cases an 
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independent salinity measurement should be obtained to allow accurate cal- 
culation of temperature in conjunction with the Brillouin speed of sound 
measurement. 

B. Methods of Measuring Salinity, 

1. Mir.rnwave Method (surface only). 

When microwaves are reflected from an objecttthe electrical properties 
of that object strongly influence the reflectivity. Microwave sensors have 
been proposed for use in SEASA’’ and were used in Sky lab for oceanograpnic 
purposes. Measurements of t.he electrical conductivity of the surface layers 
of the sea would allow the calculation of the sea surface salinity. Such 
systems have been flown in aircraft and give values accurate to one part m 
1000. We are, however, more interested in determining, salinity in depth. 

2. Raman Scattering 

The Raman spectrum of a molecule can be understood as scattering involving 
a quantum jump between two energy levels of the molecule. Since the energy of 
the scattered photon is equal to the energy of the incident photon minus the 
energy difference between the final and initial states of the molecule, the 
wavenumber shift J is obtained directly by 

J = Ei - £f 

with all energy levels in cm“^. 

Raman Scattering produces spectral bands of lower and higher energy 
referred to respectively as Stokes and anti -Stokes bands. The Stokes band 
is the most intense and is the form of Raman scattering in which wa are chiefly 
interested. The Raman effect arising from the vibration of water molecules 

is depicted in Fig. IV. 2. 

a. Raman Depolarization of H 2 O 

It turns out that both the wavelength and the polarization of the water 
Raman spectrum are temperature and salinity dependent. A recent study by 
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Chang and YoungS has shown that because seawater may absorb certain colors 
in unpredictable quantities, the wavelength shifts of Raman spectra are not 
useful in measuring temperature and salinity. However, the study further 
shows that, if the exciting light is circularly polarized, the Raman scat- 
tered light will show depolarization which is dependent on salinity and 
temperature. The circular depolarization, cannot yield the salinity or 
the temperature separately. Chang and Young proposed to treat the ocean as 
a body of constant salinity and thereby obtain the temperature profile from 
one measurement. If this method is instead combined with the Brillouin 
method, we will have two independent aquations available from which the two 
unknown variables, temperature and salinity, may be simultaneously calculated. 
The discussion of this system of simultaneous equations can be found in our 

previous report.^ 

One possible instrumentation system to utilize the Chang and Young 
results consists of a blue green laser to obtain maximum penetration of the 
sea. The laser light is circularly polarized by a Brewster window and a 

quarter wave plate. 

The Raman backscattered light from the sea is passed through an inter- 
ference filter to limit background radiation noise. The light is then passed 
through another quarter wave plate, a beam splitter and two linear polariza- 
tion filters so that the right and left hand polarization states may be 
simultaneously observed. A photomultiplier and amplifier must be provided 
for each channel. The schematic for the system is given in Fig. IV. 3. 

b) Determination of Salinity from Raman Spectroscopy of Sulfate Ions 

in the Sea. 

A ™«urenient of the ooncencratitn of sulfate lohs ih seawater allows 
the direct calculation of the salinity. >5. i’ The concentration of the 
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ions may be measured by observing the ratio in amplitude between the Raman 
scattered light and the light incident in the sample volume. 

The instrumentation required to measure the sulfate Raman effect are 
an appropriate interference filter, a photomultiplier, and an amplifier. 
Information on the intensity of light incident on a sample volume may be 
obtained from the Brillouin measurements. 

Other investigators have studied the Raman spectrum of sulfate in 
the sea, and careful attention will be paid to their results^®*^^. The 
combination of the salinity thus obtained with the independently measured 
speed of sound and depth information would allow the calculation of 
the temperature of the sea with depth. 

c) Discussion. 

The Argon Ion laser on loan to the University of Miami Laboratory for 
Optics and Astrophysics from the Navy (via NASA) will be useful in examining 
the relative merits of the different Raman spectroscopic methods. Experi- 
ments will be performed to determine which method provides the best means of 
obtaining an independent measurement of salinity. H 2 O Raman backscatter 
depolarization and SO^ Raman backscatter intensity will be measured for 
solutions of salts and typical seawater in the laboratory. 

These experiments will serve to verify the results of other researchers 
and provide experience in the techniques for future application in a com- 
bined Bri 1 louin-Raman instrument. 
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One of the reasons we are so confident of the practicality of the 
Brillouin spectrum method of temperature sensing is that the theoretical 
basis of the phenomenon is well established and that computation of the ef 


feet is possible. Brillouin spectroscopy is one of the rare fields of physics 
where theoretical work and computations normally agree with the experiment to 
within the experimental error. This agreement can be shown for the intensity 
of the Brillouin peaks as a function of the thermodynamical quantities, width of 
the peaks due to sound absorption, etc., and Fig. V.l, taken from the work of 
C.L. O'Connor and J.P. Schluff ,“°shows the agreement for the measurement we 
plan, the speed of sound in water. In this figure, the speeds of sound computed 
from laboratory observations of the Brillouin spectrum are plotted as open 
circles as a function of temperature and compared to the standard ultrasonic 



sound 3 p 66 d , indiccitsd by th6 bldck curv6* 

We use the Brillouin theory extensively in our computer simulations because 


it allows us to compute a large range of problems. Since it is the fitness 
of theory which is at the basis of our confidence in the workability of the 
design, we thought it worthwhile to study the derivation of the equations, to 
determine the assumptions used, and thus reveal any limits of the validity of 
the models chosen. The flow diagram of the derivation, seen in Fig. V.2, is 
a proof that the starting equations are basic. These equations, .iaxwell s 
equations (V.2-7) and the differential fluid equations (7.3-10), are applied 
to the solution of a specific problem: the electromagnetic waves driven by 

the polarization induced by the electric field of a strong laser beam (7.1) in an 
isotropic medium whose index of refraction is written in a T p s system. We 
emoloy various mathematical methods: Green's function (7.11), Fourier transform 


(V.l 2) and Wiener spectrum (7.13) with a rather de 


licate substitution 
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Polarization Induced by 
^of Laser Beam: 


Maxwell's Equations and Constitutive 
Equations in Polarization 

'.I 

Wave Equation Driven by Polarization 

• 

t Function^""^ 


Far Field Simplification 


■Isotropy i = a 
& 

a * n^ in 


n(T,p,s) 



Differential Fluid Equations for T,p,s 


Fourier Transform 






Wiener Spectrum 


Entropy as statistical weight of energy levels 
is thermodynamical measure of fluctuations dT.dp.ds 


Isotropic scattering coefficients in terms of macroscopic 
quantities n(Tps.x) (a^. r.. 01 

o(Tps) . . 
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POYNTING'S THEOREM IN GAUSSIAN UNITS 
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DIFFERENTIAL FLUID EQUATIONS 
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L Kx) = S(x) 

Define Green's Function G(x,x') so that 

L GCx.x') = <sCx-x'] . 

And 

i»(x) * 1 SCx' ) G(x,x' ) dx' 
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GREEN"S FUNCTION FOR WAVE EQUATION 
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FOURIER TRANSFORM 
V.12 
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Random function fUl 
of Fourier Transform FCj) 
has a Wiener Spectrut.. defined by 
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The Wiener Spectrum is sometimes called 
Power Spectrum and is used to obtain the 
average effect of random fluctuations. 
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WIENER SPECTRUM 
V.13 
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of thermodynamical for statistical quantities. These lengthy computations 
yield the characteristics of the isotopic scattering by the equations of 
Table IV. 1., For each kind of isotopic scattering, Rayleigh, Brillouin and 
salinity, the characteristics appear as these quantities: ? scattering 

coefficient, r width and v shift. These are computed in terms of the 
macroscopic quantities in the water: 

n(T,p,s,x) » index of refraction 

d(T,p,s) » density 

C » specific heat at constant pressure 

P 

V * activity coefficient of the salt 

< » thermal conductivity 

n = viscosity 

0 = diffusion coefficient of the salt 

The computation requires some first partial derivatives of the index of 

refraction and density. To this exactly computed isotropic scattering, we 
add in constant proportion the anisotropic scattering and the assumed tur- 
bidity t with the help of the turbidity normalization function. Thus, for 
each T,p,s,t we can compute the intensity scattered from the 
(15 quantities, 10 non-zero) as shown in Table V.2. 

As explained in our January 1975 Final Report, CR-139184, neither the 
salinity nor the anisotropic scattering are of appreciable intensity and, 
therefore, can be neglected. 

An example of the modeling of the spectrum of the scattered light 
return from a laser beam is shown for 20 and 40 C in Figs. V.3 and V.4. 

The separations are, respectively, of .3036 and .3114 cm *. The Srillouin 
peaks have a finite width, but it is much smaller than the separation so 
that it should not interfere with the measurement. 
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Table V.l ISOTROPIC SCATTERING 
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To the isotropic scattering 0] 

[jg p°Cai. TgCal. vgCaU 

Cag p°C9l, TgCai. 0] 

We add the anisotropic scattering 

(unimportant since r^») 5 ^ 

d 

and assume turbidity t sCah 0 Q] 

Thus for each (T.p.s.t) of the water 

and of laser with pulse intensity LG-'l 


we define Shape Function 
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and obtain the intensity scattered by a length L of water 
in a solid angle n as 


I(v) = I.(x>) * S(v.e)-^t . 


TABLE V.2 
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LASER PARAMETERS 

100 mJ/pulse 

variable 5300 A * 4000 A 

0.05 cm“^ 

4 nsec 


Table VI .1 


Power 

Wavelength 

Line Width 
Pulse Length 
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1) Yag-Nd pulsed laser frequency doubled to 5300 A: 

This laser appeared at the start of the aesign phase to be the inost in- 
teresting system: There are no limitations in power, pulse width or line 

width with suitable accessories, and powerful lasers of this kind have been 
routinely flown in airplanes. However, it turns out that the laser wavelength 
5300 A is not the best for maximum penetration in the sea, and that the acces- 
sories for reducing the pulse length and spectral width are expensive in funds 
and energy. Its use would then be reserved for a possible expensive large 
system. 

2) Oye Laser, excited by Laser: 

Oye lasers have many advantages, such as complete tunability in the wave- 
length region of interest, a natural pulse length of about 5 nanosec, and the 
additional availability of the powerful UV N 2 laser pulse itself for Raman 
salinity measurements. They however require a double Fabry-Perot in the laser 
cavity for spectral width reduction and are at present limitec in power by 
dye destruction.. Since they have efficient energy conversion, they may be 
used for the airborne system, if the technology of .more powerful large bore 
lasers with narrow spectral width continues to make progress. 

3) Argon Ion Laser 

These provide easi'y sufficient continuous power in a series of very 
narrow spectral lines in the region of interest 50QG - iSCC A with a well proven 
technology. Hcwever, the need for modulation in the monostatic systems will tend 
to slow down the measurements and reduce the power availaole. This laser sys- 
tem clearly appears to be the least expensive and the- easiest to use for the 
first phase of the experimental program. 

3. Optical System 

The optical system, as shown ’"n Fig. 71.1, may be kept rather simple: A 

method to bring the laser oeam to the water to be analyzed, and a large lens 
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for collecting the scattered light. Since all considerations indicate 

definite advantages of keeping the laser beam as small as possible, we nave 
only to steer the beam by a mirror, and the only attenuation is aue to the 
absorption h(x) of the underwater path Y (the small reflection at the surface 
should be negligible and will be included -in the- sea surface factor). 

The scattered light is then collected by the largest possible lens (of area A) 
at a height H above the water. Taking into account the refraction at the sea 
surface, the solid angle collected is ^ where n is the index of 

refraction of the water. The power scat‘Ured in this solid angle is then the 
laser power multiolied by the scattering coefficient, the solid angle, and the 
length of the scattering medium observed. This length L is found from the 
geometry in a bistatic system and from the integration time of the detector 


,t, in the monostatic system ^ '' 

attenuated by the absorption of the underwater path by a rac.or e and 

by losses at the sea surface described by a factor k, a function of the dimension 
of the re- . beam. In conclusion, the light intensity pa-ssing through the 

area A of the lens can be computed as a function of wave number. 

cA . 

I(^;) = k {Iq(v) * o(.,a):- 2r,(Y + nH)'- 

and the instrument has to make the best use of this ocserved lignt. 

C. Dispersive System 

This collected light has to be analyzed for the spectral information tnat 
will give the speed of sound and the turbidity. This is done by distributing 
the photons in several channels; to obtain two quanti-.iss . Oi 

turbidity) by ratio of intensity in the channels requires tnree cnanre.s. 

After discarding seven'll interesting schemes or diiOers./e a/s.e. s. 


rhis scattered intensity is 

... . -h(A)Y 
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we are now considering only two options: Fabry-Perot and Mach Zehnder inter- 

ferometers. We will probably use the Fabry-Perot system in tne first stage 
of the experimental work because a suitable Fabry-Perot interferometer is 
available at our laboratory and we have years of experience in its use, but 
because the Mach Zehnder system is not as sensitive to the effect of the sea 
surface and is more stable, we plan to use it for the airborne system. 

Fig. VI. 2 illustrates the principle of the Fabry-Perot: Each wavelength 

A is transmitted in a ring of angle e. A multimirror inclined at an angle j 
can then direct various wavelength ranges towards photomultipliers at angles 
2o (Fig. 71.3). The quality of tne Fabry-Perot is best expressed by the 
"finesse" F, defined as the inverse of the unavoidable spectral spreading 
compared to the free spectral range FSR. The effect of the Fabry-Perot can 
then be written as the convolution of the spectra with the spread function 


S(., :i 

FSR2 + F^v + k FSR): 


and the multimirror 


selects a suitable wavelength range in this convoluted 


spectra. The range of the three channels are shown in Fig. VI. 4 and VI. 5 
compared to the convoluted Srillouin spectr'a at 20 and 40°C. The signals 
are given for a shipborne system looking at a depth of 20 meters with a 0^1 
joule laser at 4500 A in units of 10' photons. The schematic of the Fabry- 
■^erot system is given in Fig. 71.6. 

The Mach Zehnder system is based on another interfe'^cmeter shown in 
Fig. VI. 7. This instrument is often thought as an unfolded Michel son 
i nterferometer yielding the same channeled spectrum, but in reality the un- 
folding allows the transmission of two compl ementary channeled spectra, 'he 
advantage of a Mach Zehnder is a consequence of tne fact that no protons are 
rejected: This not only results in a great sensitivity, but it also frees 

the design of the necessity of insuring uniform illuminat'on o* tne channels. 
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The amount of light accepted by the interferometer is easily increased 
by field widening as in Fig. VI. 8. To obtain the necessary three channels, 
two Mach-Zehnder interferometer are used in succession: A first one is used to 

separata the Tyndall peak, and the second one for analyzing the remaining 
light to get the Srillouin separation. The transmission in these channels 

will be of the form 

T(v) " a[l - 3COS (1^(.v+Vq)][ 1 ± YCOSa2(u+'Jg)] . 

A schematic of the complete system is shown in Fig. VI. 9. and the expected 
transmission of the Brillouin channels in Figs. VI. 10 and VI. 11. 
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0. Collection and Recording 

The systems used to detect and record the photons in eacn cnannei 
have to be quite dirfarent in a bistatic, monostatic oulse, and monostatic 
modulated Srillouin remote sensor. The bistatic configuration is certainly 
the simplest because continuous or chopped light has to be detected; We will 
simply use good sensitive photomultipliers with quantum efficiency in the blue 
of about 20 to 30%. With these detectors we don't expect to have to use a photon 

counting mode even though this option is available if more sensitivity is 
desired. The voltage signal of the three channel will then be read sequentially 
by the NOVA 300 computer. The material for the analog data acquisition system 
is already on hand and we are wor!<ing on the modifications to the comouter. The 
data accumulated in the computer during a complete run will then oe dumped on 
teletyoe or magnetic tape for future reference. 

ihe detection system is a dirricult and expensive part of tne monostatic 
pulsed system. The extremely fast varying signal must be measured precisely 
with a time resolution of a few nanoseconds; just maintaining a pulse shape 
in the electronics at these speeds is a complex undertaking. The best solution 
is to reduce the electronics to a minimum and thus use fast photcmultiol iers with 
large gains (tnese may have long delay times but are designed to avoid the spread 
and jitter in this delay) and then, after minimum amplification, we digitize 
and slew down the pulse by the use of a Textronic trans'ent dig; tizer, wnicn 
can digitize signals of mVol: range at nanosecond speed. This however represents 
a large expense (about 33,-iQO per transient digitizer), and we will need 3 digitize 
(one per cha.nnei) and a soeciai control computer. For tasting purpesas, a c.neacer 
systam without continuous depth information could procaoly be rigged uo ^rqm snort 
aperture sa.mplas and nolds. This may be a met.ncc 'or measuring tne :ua’’ty n- t:.e 
onotcmultiol iers without the exoense of a ful] recording systam. 
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E. Data Processing 

For reasons of convenience, we will record the data in a minicomputer, 
which allows great flexibility in subsequent processing of the data. However, 
a relatively simple linear algorithm such as shown below will be sufficient for 

the prototype; 

The intensity in the three channels described in Section C of this chapter 
are assigned the notation I3 in the following way: 

I. is the channel that increases most with an increase of 


Tyndall scattering (central channel). 


I and U are the channels that resoectively decrease and 

2 j 

and increase with an increase in the separation of tne 
Brillouin peaks. 

Then in a large number o'‘ experiments, either simulated or real, the temper- 
ature of equality of I2 and I3 is determined, and a regression on data 
at this temperature but at different turbidities gives the coefficients and 
of Tyndall rejection by requiring that the quantities 

I4 = I2 " I3 ■ Vl 

h ^ ‘3 ■ '2 ^ ''^^i 

became independent in first order of the turbidity. 

It then becomes easy, for this case of assumed salinity, to rind tne 
coefficients for the linear expressions for temperature and turbidity by 
further regressions: 

t = I1/I4 

T = T^ * C3 I5/I4 . 

This linear processing was used in the simulations for the design 0. -i ■ - 
ferent configurations, so that the noise could be obtained by a method similar 
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to the square root of the signal. We also discovered at that time the 
remarkable precision of such a simple computation; no noticeable error 
with turbidity from .01 to 100 and drifts of less than PC for range of 
more than 20°C. 

An easy method for avoiding the error due to even this drift is to 
make a table of the needed corrections as a function of the raw resuPs, 
such a correction table may already be stored in the computer memory for 
the bistatic system so that the direct result of an experiment will be 
given immediately by the computer as correct temoerature. 



VII. Feasibility Criteria 
A. '^pahnrne 3i Static System 

The feasibility of the seaborne bistatic system is indicated in part 
by numerous laboratory experiments to measure tbe speed of sound from the 
Srillouin spectrum. Typically, a laser beam is directed into t.be liquid, 
and a Fabry Perot is used to analyze the spectrum. This corresponds almost 
exactly to the prooosed bistatic system except that we will encounter field 
conditions and impure water. 

Since we demonstrated the overall feasibility ot the Srillouin tempera 
ture remote sensin,. -e feel that the intensity of the return and freedom 
from noise will be quite sufficient in the bistatic system. Comparinq with 
the detailed computation on the saaborne monostatic Fabry-Perot system discussed 
in this Chapter, using a 1 second integration time with a 650 m« laser which 
corresponds to 650 mj instead of only 100, and assuming that the bistatic 
system could easily look at the same scattering lengths (about 1 meter), 
we conclude that measurements up to a depth or 20 meters should be possible 
with the limited instrumehtation we are contemplating. Our confidence results 
in good part from our previous experience with the rabry-Perot ihterrerometer 
and nlgn sensitivity photomultipliers. «e have developed techniques for 
modeling and assessing these systems, ae know that if extreme sensitivity is 
needed we can use photon counting; and with the resulting data in a computer, 
many processings and time integration schemes may easily be tested. 

E/en without searching at this time for the exact capabilities of a 
working bistatic system, we will need the experience that such a simple system 
with flexible design only can bring. It not only is the si...ples, .nd os. 
direct demonstration of the principle of the Srillouin temperature remote 
sensing, but also offer a similarly simole and 1‘ract method of c-eckinc tne 
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whole of the theory we have developed. The verification can start with a 
situation vsry close to ths normal laboratory Brillouin sxper’msnt, and than 
be varied continuously up to a practical remote sensing situation in 
biologically active water with normal turbidity at a depth around 10 meters. 

B. Seaborne Monostatic System 

This system was completely modeled in a computer program on the Nova 
800. The expressions relevant to the simple system: laser, scattering, 

collection lens and Fabry-Perot were taken from chapter V and 71 of this 
report and programmed in Sasic. A printout of the program is given in 
Appendix 3. The results appear as the typical computer run shown in 
Tables VII. 1 and 2. After a series of questions and answers about the kind 
of instrument to be modeled, the computer determines the best coefficients 
of the processing on the data, and thereafter, for each temperature, salinity, 
turbidity and depth, the program computes the expected signal of each channel 
in number of photons. These channels responses are orocessed to obtain the 
temperature and turbidity, and their root mean square variance. 

The most interesting results of this study are expressed in Fig. 'VII. 1: 
Temperature precision in degree centigrade (V) compared to laser wave- 
length for differ'ent depth, and Fig. VII. 2: Turbidity precision m per 

cents (Ic) compared to laser wavelengths for dirferent deptns. .re 
program was also used to give numerical answers to t'wo imoortant ouestions 
on the system: The influence of the turbidity on tne measurement, ar.d the 

effect of an error in salinity. The most noticeable e*fect of turbidity was 
on the precision in temperature, and is shewn in Fig. VI 1. 3, as a runction or 
the finesse of the Fabry-Perot. It is however a very encouraging sign that 
this effect begins to appear in relatively dirty water; A turbidity or 13 
corresponds to half transmission lengtn of only 52 cm. .^s seen in -ig. V.I.d, 
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the effect of an error of one part per thousand in salinity corresponds to a 
drift of 0.6 ‘’C. 

These computations give us full confidence in the practicability of a 
seaborne monostatic system, but there is no doubt that utmost care will have 
to be exercised in the design of the detection and recording system. Fortunately, 
fast light pulse detection has been the subject of many tf-chnological studies, 
and our requirements for a monostatic system are within commercial specifications. 

C. The Airborne Monostatic 

The feasibility of the airborne system was demonstrated in tne NASA reoort*- 
CR-1 39184, January 1975, and the computations made so far have only increased 
our confidence in those computations. Since then, we have improved tne design 
of the airborne system by the proposed use of a douuie Mach-Zehnde'' interferometer , 
thereby increasing the luminosity by a factor of more than three and avoiding 
troublesome sea surface ef‘*'ects. 


0. Extraneous Light Moise 

In a chapter on feasibility, it is adequate to discuss the problem of tne 
possibility of error arising from extraneous signals in the photomu i tipi ier . 

All our computations indicate that with suitable interference filters, there 
is no extraneous light p'-oblem in a monostatic 3rillouin temperature remote 
sensor because the scattered lignt is distributed in an extremely small soectrai 
and temporal resolut’on. Actually, even assuming .nat L..e sea sanu. 

all the -eflected sunlight possiole into the instrument, a Fabry-Perot working 
at spectral resolution 0.02 A and temporal resol uticn of 5 ns, receives only 
extraneous light o*' 1,300 photons to bn compared with tne exoected 10 

of signal . 


crotons 


VIII. Bistatic System 


The next phase wiil be to prove directly that the Brillcuin scattering 
method can be used in a practical way to measure temperatures outside the 
laboratory. A laser will be mounted on a barge or dock and its beam will 
probe the water to a distance depending on the water clarity at the site 
avai lable. 


In detail the bistatic system will consist of a platform which will 
initially be a dock, and ultimately, if possible, a ship, upon which is 
mounted the technical equipment. The instrument is shown in Figure Viil.l 
and consists of a laser system and a receiving system. 

The laser system will consist of an RCA-LD 2101 Argon Ion Laser 
equal) complete with power supply. There will also be a mirror system 
with a flat under-water window to guide the laser beam into the water. 
Adjustments for aiming the laser beam will be orovided. 

The receiving system will include a second mirror system with lens and 
under-water window to collect the scattered light from the sea. ^e have 
available for comparison a 9"x24" plastic Fresnel lens, a 6 achromaw, and 
a 10“ apochromat). An optical train must be built to carry the light to 
the Fabry-Perot i nterferometer . The interterometer itself, which is already 
in the laboratory, must be adjusted and mounted. A triole mirror C'Fafnir" 
in Fig. VI. 6) divides the spectrum into three model 95536 ohoto- 

mul tipi iers , as shown a’so in Fig. VI. 5. These photomul ti pi lers as weii 
as their power supplies are already in our laboratory, but trey mus, ..e 
mounted in place. 

The signals froti; the pnototubes will be ampli-ied by tnree model 1225 
Teledyne =hilbric'< ~l~ ooerational ampli*-ers, and tneir cutout *ed tnrcuc.n 
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3 charnels of an 3-channel data acquisition "Hyorid Systems model 0AS400" 
into a Nova 300 computer. The amplifiers, 3 channel system, and computer 
are already in our laboratory. 

The raw data will be available in a magnetic tape dump, but the 12 
kiloword memory of the computer will allow quite complex processing in 
real time and the resulting temperatures and turbidities will be available 
on an ASR33 Teletype, also in our laboratory. 
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IX. Implementation Plan. 

A. Bistatic Experiment . 

Although the ultimate airborne system will be monostatic, with essentially 
parallel laser and scattered-light beams, the next phase of this work will be 
bistatic (see flow chart. Fig. IX. 1), with separated sendiV.g and receiving 
locations. The reasons for using this intermediate step are twofold: 

1) The budget for the period March 1 , 1976 - February 28, 1977 was cut 
from the requested level of about $100,000 to $50,000. This meant that the 
amount to be spent had to be greatly reduced, and that the experiment will be 
based on materials already on hand or borrowed, or on surplus eouipment. 
Accordingly, we plan to use a continuous argon-ion laser borrowed from the 
Navy, and electronics already in our Laboratory. This necessitates a bistatic 
system. 

2) It is probable that even without the budget cut, we would have initially 
used a bistatic method. This is because it is simpler in execution than the 
monostatic system. Temperature and turbidity information at depth have never 
been obtained by laser teams in the open, and it is best to begin with the 
simplest possible configuration. The vital thing is the demonstration that 
practical laser beams in the sea can indeed generate the Brillouin, Raman, 

and Tyndall scattering in intensities sufficient to be useful. 

This can be demonstrated as well with the bistatic experiment planned 
for the next phase as with the considerably more complicated monostatic method. 
The development of the systems through bistatic, seaborne monostatic, and 
airborne is shown schematically in Figure IX. 2. 

3. How the Monostatic Method 2eve1cc.> f’^om the 3i static . 

For airborne measurements the effect of the sea surface can be min'miced 
by utilizing the same optical pa tn for dcwngoing and uocoming oeams. This means 
t.hat tim.e of flight, ratner than the angle, is used for deotn ’n'ormation. Sine 
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FIGURE IX. 





one location is used for projecting the downgoing beam and receiving the 
upcoming beam, this is called a monostatic system. 

The physical principle, however, with the exception of the depth 
measurement, is exactly the same with the monostatic as with the bistatic 
method. Srillouin scattering is used for obtaining the speed of sound in 
exactly the same way. A combination of Tyndall and Srillou’n scatuering is 
used to obtain the turbidity, and finally the Raman scattering is used in 
the same way to measure the salinity. Only the depth^measurements differ, 


angles being used with the bistatic system and time gates being used wi-h .he 
monostatic method. Since the principle of the experiment is almost identic?:, 
the monostatic seaborne program follows logically from the bistatic program 
and prepares directly for the airborne work to follow. 

There will, however, be equipment changes, as explained above. Since 

depth informiation is to be obtained by tne time delay between tne dcwngoing 
laser and returning scattered light, a pulsed laser or a modulator for tne 
laser beam will be necessary. A fast electronic receiving and recording 
system will also be necessary. All these comoonents are now available *rom 
commercial establishments. Pulsed lasers and modulators have long ..een in 
use, (in fact a pulsed ruby was the first operating laser made) and recently 
Tectronix Incorporated has brought out a suitable detectcr-recordi ng iys.em. 


C. How the Ai’-oorne S/stem Develops ~rcm tne ‘^cnostatic Seaocrne oystem . 

The basic requirements of an airborne system are adecuate pcvver, ccmspacuness , 
lightness, and reliability. In principle, tne system is identical with the sea- 
bo’-ne monostatic 3>stam. Therefore, our experience «ith the seaborne system 
should provide c-ucial information on the choice of tcnocnents *or tne amocrne 
system. Because of tne weignt and compactness '■ecu' rements *n an a rcra ,, 
we probably will have to custom-design tne power suoply sections or tne 
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apparatus. The laser itself, together with the receiving and recording 
electronics, will probably be identical with or only sligntly modified 

from the sea-borne system. 

0. Design Differences Between Stages 

The design changes through each stage are shown in Figure IX. 3 , and 
the following discussion refers to that figure. 

1) The basic method (except depth measurement) is the same throughout. 

2) The dispersive element for tlis first and second stages is planned to 
be the Fabry-Ferot. This is because a suitable Fabry-Perot interferometer 
is already in operation in our Laboratory and has been used successfully 
in preliminary experiments on the Brillouin effect itself. As explained 
elsewhere in this report, however, a tandem Mach-Zehnder is the ideal instru- 
ment for the monostatic airborne phase, and it will, therefore, be used at that 
point in the program. 

3) The optical system will be somewhat different in each phase. In the 
bistatic configuration, the projecting and collecting optical systems are 
separated by the baseline distance. There must also be • means for changing 
the angle between the two beams. The two monostatic configurations will be 
almost identical with col inear beams, but because of the more stringent 
space and weight requirements in the airborne phase, there will probably be 
some redesign called for at that point. 

4) The bi static detector will be a standard high-sensitive photomultiplier, 
such as we have already in our Laboratory, since the bi static signals will be 
low-freguency in nature. For the two monostatic phases of the work, a high 
speed photomultiplier and associated electronics will be used. 

5) The power supply for the bistatic phase will be the unit already in our 
Laboratory, whi ch was delivered with the RCA argon-ion laser. In the two 
later stages, power supplies for the pulsed laser will be necessary. Although 
the same laser is planned for both of the monostatic phases, the space and 
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weight requirements of the airplane will probably necessitate a special 
version of the pulsed laser power supply. 

6) The recording system for the bi static phase is already in the 
Laboratory. Since time resolution is not used to measure the depth, the 
system is essentially d.c., and the problem is easily solved. In the two 
monostatic cases, a very fast recording system is necessary. As mentioned 
above, at least one commercially available unit CTectronix) is available 
and suitable, and it will operate both for the seaborne and airborne phases. 

7) For the bi static phase, our plan is to use the RCA two-watt argon ion 
laser already in our Laboratory. As discussed elsewhere, a pulsed or modulated 
laser will be necessary for the two monostatic phases. We plan to use the 
same laser for both the seaborne and airborne monostatic phases to cut down 
expense. 

£, Cost and Time Estimate of tach Phase . 

Because of the'continuous development of the state of the art, especially 
in the areas of powerful pulsed laser and data handling equipment, our 
estimates for Experimental Phases II and III (nonostatic, seaborne and air- 
borne) are considerably less certain than for the first Experimental Stage. 
However, it is possible to make rough estimates of the cost of principal 
subsystems. These are listed in Table IX. 1 and Table IX. 2 below. 

The time for each phase is now estimated to be one year. See Tables 
IX. 3 and IX. 4. 
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TABLE IX.l 

TABLE OF VALUE OF COMPONENTS 
BISTATIC SYSTEM 


Approximate 


Item 

Source 

Value 

Laser System 

Loan from NASA 

$16,000 

Window and Lens System 

To be purchased and 
constructed 

1 ,000 

Fabry Perot 100 mm 

Laboratory 

6,000 

3 EMI 9558B Photomultipliers 
and Power Supplies 

Laboratory 

3,000 

Operational Amplifier 

Laboratory 

1 ,000 


Teledyne Model 1025 
and Power supplies 


Data Acquisition System 
Hybrid DAS400 

Laboratory 

500 

Nova 800 Minicomputer 
with 12 K Memory 

Laboratory 

8,000 

Teletype ASR 33 

Laboratory 

2,000 

Magnetic Tape Dump 

Laboratory 

2,000 

Trailer for Site 

To be converted from 
surplus 

'2,000 

Power cables and 

To be purchased and built 

1 ,500 


connection for si te 


Total 43,000 

of which 38,500 

is already on hand 


TABLE IX. 2 


COST BREAKDOWN (APPROXIMATE) OF ENTIRE PROJECT. 


I. 


Experimental Phase I March 1 976-February 1977 

Salaries, Benefits and Indirect Costs 
Supplies, Equipment, Reproduction, 

Computing and Related Items 
Travel and Site Expenses 


36,000 

7,800 

6,200 


TOTAL 

II . Experimental Phase II May 1977 -April 1978 

Salaries, Benefits and Indirect Costs 
Supplies, Equipment, Reproduction, 

Computing and Related Items 
Travel and Site Expenses 


TOTAL 

III. Experimental Phase III May 1978-April 1979 

Salaries, Benefits and Indirect Costs 
Supplies, Equipment, Reproduction, 

Computing and Related Items 
Travel and Site Expenses 


TOTAL 


$50,000 


60,000 

30.000 

10.000 

$ 100,000 


80,000 

60,000 

10,000 


$ 150,000 


TOTAL FOR THE THREE PHASES 


$300,000 


mental Stage 



table IX. 3 
















eriinental 
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X. Summary and Conclusion 

A, Hnw the Parameters Are to Be MgASMTSd- 

As detailed in the Report above, the oceanographic parameters to be 
measured are sound, velocity, turbidity, temperature and salinity, all as a 
function of depth. The measurement principles are shown in Table X.l for 
each phase of the Experimental Program. Mote that as far as the principles 
are concerned, there is no difference between the phases except that the 
depth measurement in Experimental Phases II and III is by time of flight while 
in Phase I it is by trigonometry. 

Of course as we move from bi static to monostatic seaborne to airborne, 
the requirements on equipment that must be used become more and more stringent 
in weight, space, and reliability. This means that the components 
are continually increasing in sophistication as we proceed through the 
phases . 

B. Criteria for Determining Feasibility and Acceptabi.lity. 

In any program where something is to be measured, the accuracy of the 
resulting measurement is a primary goal. Accuracy of measurement of the various 
parameters in the program is therefore a criterion for feasibility. In the 
feasibility study, we found that according to our calculations, the accuracy 
expected from these methods depends on the length of time of the observations, 
the depth of the observations, and the clarit; of the water. The accuracy 
also depends on the design of the equipment to be used, and we concluded that 
practicable equipment could be assembled which would give useful information. 

In order to prove feasibility we should show that the laser probe method 
is superior to other methods for determining the parameters in shallow seas. 
This superiority is expected to show itself in the much greater speed with 
which the observations can be accomplished, as described in the Introduction 
to this Report. If temperature measurements at depth can be made by laser 
scattering from a ship or airplane, a very large improvement in information 
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Table X.l 
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gathering capability will have been provided. The accuracy at a single point 
cannot match the thermometer method, since thermometers can measure to a 
small fraction of a degree, but in the study of the ocean, it is the measure- 
ment of temperature gradients that is important, and here the laser method 
can disclose infonnation that the point-by-point thermometric methods would 
not obtain, either by providing too coarse a grid or simply by taking too 
long for each observation. Therefore the criterion for acceptance is whether 
or not useful data can be gathered by our method that would not have been 
obtained by the classical methods. 

c. Proof that Raman and 3rillouin Sca ttering Can Be Measured Remote I jj. 

Of course, the best proof is the actual carrying out of the program. 

This we hope and plan to do. In the meantime we must base our judgments on 
the best information available. These are listed below. 

1 . Theory 

The Theory of Srillouin Scattering is ore of the best examples of the 
application of the principles of physics. It is described in detail in Chapter 
II of the Detailed Technical Report of our Feasibility Study, NASA CR-139184 
dated January 1975, and is developed in Chapter V of this report. This theory 
has been checked experimentally in many laboratories over a period of almost 
forty years and has been found correct within experimental error. We have 
made laboratory studies of our own described in the foregoing report 
which agree with the theory in every way. The theoretical material relat.ng 
to the optical part of the plan is also beyond question. Our own experience 
with similar methods (see the publication lists in the Aopendix to this 
Report) have provided very high confidence in our predictions. 

As we have stated in the body of this report, the calculations we have made 
borne out by our experiments in the Laboratory, show that field measurements 
of the oceanic parameters can be made. 
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2. Experiments by Others 

Many of the elements that make up our Program have been performed in the 
field by other groups, for various purposes, greatly adding to our confidence 
that all of these measurements can be accomplished together. 

a) Depth measurements have been made routinely from aircraft by laser. 

Just as in our system, a low-flying airplane was used as a platform for a 
pulsed laser firing vertically downward. Light scattered from the bottom 
was received back in the airplane and the time difference used to measure 
depths. Accuracies of at least .5 meter were considered routine, and most 
important, the state of the sea (providing that the conditions were reasonable) 
seemed to have minimal effect on the accuracy of the observations or the 
intensity of the back-scattered light. 

b) Raman spectra have been used to measure impurities (.oil spilis) on 
the sea surface by a laser probe from an airplane.'^ Here again, the fact 
that Raman scattering, which is of the same order of intenstty as the Bril- 
louin spectrum wt plan to use, can be observed from an aircraft is very 
encouraging indeed. 

c) Brillouin scattering has been used routinely in the laboratory for 
the measurement of sound speed in liquids. Hare again, the fact uhat it has 
been found desirable to use this method is an indication of its practicability 
that argues well for the success of the program. 

3. Conclusion 

Thus, our conclusion is based both on theory and experiment, and nothing 
we have learned has decreased our confidence in the program. In fact, as each 
new result has been obtained, our expectation of success has increased. 

0 . A Brief Description of the Following Phases of the Program . 

1. The next phase will be to prove directly that the Brillouin scattering 
method can be used in a practical way to measure temperatures outside the 
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laboratory. A laser will be mounted on a barge or dock and its beam will 
probe the water to a distance depending on the water clarity at the site. 

Because of fund limitations, this may be in the Miami area initially. 
Scattered light will be received at a nearby location and the Brillouin 
spectrum analyzed aS function of distance. This is the bistatic experiment 
described in detail in the body of the Report. Success here will prove that 

the method is practicable in the field. 

2. The following phase will employ the time of flight method of deter- 
mining distance, which is the method that will be used in the final, 
airborne phase. The results of the preceding step will be used he-^e, but 
now the laser is to be pulsed in order to provide time data. Recording 
instrumentation for this purpose must be added at this point. 

3. The final stage is airborne, as described in the body of the Report 
(see Figure 1.6). Here, the space and weight requirements of an aircraft 
will necessitate a certain amount of redesign. Also, the equipment will be 
partially automated, and special care will be concentrated on questions of 
rel iabi 1 i ty . 

4. At each stage, reports and visits will be used to insure close 
cooperation between the Laooratory and NASA. A flow plan for the entire 
project is shown in Figure X.2. 



FIGURE X.2 

FLOW CHART FOR ENTIRE PROJECT 
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PROGRAM FOR SEABORNE MONOSTATIC SYSTEM 


10 
15 
I 7 

aa 

ai 

as 

26 

27 

23 

32 

33 

35 

36 

37 
33 
39 
43 
4 I 
4a 
43 
44 - 

45 

46 
50 
55 
ST 
60 
61 
6a. 

63 

64 


PRINT ••LASEH:F0W£R>LA«B0AvVpTH^PUL3S. L£NGiiT” 

INPUT P1.»LjW7.»NI 

?rIn? «0BSEHVAT ion; height, pupil AHEA,3£A surface FACTOR,”; 
PRINT "EFF IC lENCY, ANGLE” 

INPUT H,A,HliKa,Ta 
LET H*100»H ' 

PRINT -SINGLE F.P* INSTRUMENT . G WE FSR, F INESSE,OAMP ING TIME_ 

GOSU3 "7500 
LET 31*1 
INPUT Fl,F3,Na 
® R INT 

PRINT -ASSUME SALEN ITY OF WATER- 

INPUT 3 

PRINT 

let T»as 

GOSUB 34 0 0 
G0SU3 3600 


GOSUB 

GOSU3 


3700 

3000 


PRINT •’SHIFT AT 25 C IS :'JHl 

oaiN?^-GIVE THE CHANNEL LIMITS C6 WAVENUMBERS); 
INPUT 3 C‘ 3 ),aCl l,BCa3,3C3l,aC4 3,BC51 
PRINT 

0 IM XC6 ] ,YC6 1 

let X3»6E-*-ll*L>*‘Pl*Kl*K2*A»Na/ra 

LET A3» COS a .74533E-2’*Ta) 

LET A4»3»<1 ♦A3*A3 )/4 
let A3»3»d3-*>A3*A3>/40 


6 5 

LET W3»W7i*Fl/F2 


67 

LET A3*100*‘l'3»- 

aS-T2)-^»; 

63 

GOSUa 7900 


69 

LET K3«K3»Fl/3C5 1 


70 

LET P2«l 


71 

LET AS-A5/100 


7a 

LET T 6 • 3 


75 

FOR T» 0 TO 60 


30 

GOSUB 4000 


35 

IF <12-13 )<• 0 

GOTO 1 3 

90 

NEXT T' 


95 

PRINT -ERROR; "0 

LARGER 7 

97 

GOTO 40 


98 

PRINT -ERROR; 73 


99 

GOTO 40 


100 

let la^ia^IO 


131 

IF T<3 GOTO 98 


1 05 

LET 13 -13 -12 


1 1 0 

LET 76-13 


115 

GOSUB 4000 





155 


T 


I 20 
I 25 
1 30 
t 35 
140 
145 
150 
152 

155 

156 

157 
155 

159 

160 
165 
I 70 
1 75 
1 30 
LSS 
190 
1 91 
192 
330 
305 
336 
307 
303 
3 10 
3 15 
3 1 7 
3 20 
3 22 
3 25 
327 

330 

331 

332 

333 

334 

336 

337 
333 

339 

340 
345 
347 
350 
355 
360 

365 

366 

367 


L£T C2* C 12^I3-14 )/ II 
L£T C4*'< Is ♦12-13 ll 
L£T Cl «<T6*l4 II ■ * 

L£T 76*‘ 0 

FO R T «T TO 0 STS? - • 1 
GO SUB 4000 

£p- ( X3 — X2>C4'^I1 ) ■** 0 GOTO 155 
NEXT T ■ ■ ■ 

FOR T»T TO 60 STS? *01 
GOSUa 4000 

IF < 13 - I2+C4»ll ) ** 0 GOTO 15 9 
NEXT T* ■ ■ 
let T0»T 

LET CS«C1»II/CI2-^13-C2>«^I1 ) 
let T*T3^10 ■ ■ ' 

let C3»l3*t I3»I2-C2»tl )/ C l3-12*'-4*il 


) 


"TKE‘ 

''LBT^ 


processing IS';’* 

.'♦jci ;”il/ c I2^i3^’*/*C2;’* 11 )- 


«•« 


;T3. 


;C3: 


( 11 V C 12^p“">C2J*J 11 


••GIVE ACTUAL VALUES 0? TEMP* 
r,3^T6.»P2 
«130»P2 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
INPUT 
LET P2 
PRINT 
PRINT 

GO SUB 4000 
LET 14* 12^13-C2»11 
15*13-12+C4*I1 
N4*- SaR CI2>I3+C2*C2*I1 ) 

N5» SaR ■c12V13-H:4*C4*11 ) 

T7-C1 *11/I4-C5 
T5»Cl«-'SaR CIl)/ 14 
T9*C1*N4*11/ i 14*14 ) 

TS»100* SCR' 'CT5 *T5^T9*T9 )/T7 
T9»130»<T5+T9)/T7 
T4»C3*NS/ 14 
T5*C3*IS *N'4/ < 14*14 ) 

T3* SaR CT4»T4^T5 »T5 ) 

T5»T4*75' 

T4*T3 

T3*T0*C3*I5/ 14 
P3»14*9896*<N1 ♦N2)/N 

••THE INSTRUMENT GIVES'* 

''CHANNEL'S SIGNALS OF'*/ 11 / 12/ 13 
•*’7 15; TSSPERATURE 0 F'*/T3 /'* INSTEAD 
‘*’7lTH A RMS 3ETWSEN’*7T4;'*AN0'*/T5/' 
••AND A TURBIDITY OF 
'•’U ITH A RMS ■ 3ET'7E£N 


SALlNlTY/TURSp ITY_/DEPTH'J 


LET 
LET 
LET 
LET 
LET 
LET 
LET 
LET 
LET 
LET 
LET 
LET 
LET 
LET 
LET 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 


•/T7; 
: • 


OF* 
’C’*' 

•• INSTEAD ‘Cf’ 

j i y j ^ 


A I 


s;”And 


••ano an 


IMPREC 13 ION ' In depth. OF”/ ?3/'JGM'J 


< l.< U ' ' : i 'i-i- 
Ol> 
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365 PRINT 

369 PRINT 

370 GOTO 300 
4000 L£T K4«K3= 
4033 L£T P»9S0 


^ EXP <a 7*P2)/C(P2^1 •33*H)*<?2 + 1 *33*h)) 
,663»P2+l •01325E-*-6 


4010 

GO SUB 

8400 


4015 

GO SUB 

3600 


4020 

GO SUB 

5700 


4025 

GO SUB 

3900 


4030 

GO SUB 

9000 


4035 

GO SUB 

9100 


4040 

GO SUB 

9200 


4045 

GO SUB 

9300 


4050 

GO SUB 

5000 


4055 

FOR a 

* 0 TO 

2 

4060 

LET 

xcai- 

0 


4063 

4070 

4075 

4050 

40 55 

4090 

4095 

4100 

4102 

4105 

7500 

73 05 

7510 

7515 

73 20 

75 22 

7525 

7900 

7910 

5000 

5005 

5010 

5015 

3020 

3025 

3030 

3035 

3040 

5045 

3050 

3055 

5060 

3063 

3070 

3400 

3420 


FOR U2*BC2*ai 
GO SUB 9500 
let xca3»xtai 

NEXT U2 
NEXT 2 

LET 11 *K4*Xt 0 3 
LET 12»K4*XC13 
LET I3»K4*XC21 
PRINT 
rlETURN'" 


TO SC2^‘a^l 3 STEP •301 


>.001 *UL 


C9*4 .1 3£-^T 
F9»l' 

U=.0S 
K*5 9200 
G«.02 
O«.0000l3 
RETURN 

LET A7»-.30034-S.6£”10^CU"^550 ) » CL-4S5 0 ) 

RETURN 

LET T1«T-^273.16 

H6»N#T1 »D1/ CL*L) 

R6“R6*R6*4 .36564S*! 3/ CR0*C9) 

oA M ^ »C 6»D2-0 1 »Al »T L") *N/ CL»L) 

S6aB6»B6«4 .36364Ei>l 3«T 1/ r<C9*R0 -SI -A1 -Al -TL ) »R0-C9 ) 

36-N-03/ <L»L) 

S6-36-S6-3 •49E-*'13-S/f 9 
A6-10-U-<R6 •►36+S6 )/ (6-7-U ) 

A2-4094 ►7»N« SIN C 3 »7266£ -3 -T2 )/ L 
V0-W8+K-A2-A2/ iR0'»C9) 

V 1 ■ W 3 +G -A2 -A2? "C 2 -R0 > 

V2-WS-D-A2-A2* " 

V3-WS + .31 532-Tl 
HI -1 .62924E-6-V-A2 

7 CS^T^P ) 


LET 

LET 

LET 

LET 

LET 

LET 


LET 

LET 

LET 

LET 

LET 

LET 

LET 


LET 

LET 


RETURN 
REM 

let 74-1402.3 9 


SPEED or SOUND <7 ) IN SEAWATER AS 
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t 


L£T vs ■ <S 1 1 3 9 ) *T 

LET VS »V5 -<S •50 947E-2 ) »T »T 

let VS»VS +‘<2.2lS36£-4 )»T>*T*T 
let V6»a.32952)*S»a00a) 

let V6 »C^6 -^< cl .28956E-4 )*S*S»ac.*o ) ) 

LET V7«C.l550S9)*P/ C9.8LE+6) 

LET V7»V7 + C C2.4A999E-5 )»P»P/ ^9.8L=.-*-6 

let V7«V7-C’C8.83392E-9)*P»P»P/ C9 .81C.+6 1 3 ) ) 
let V8»-aV27563E-2)*T*S»aa00'> 

let V8*VS-*'CC6.3S192E-3)^T*P/ C9.5li+o>) 

let V 8-V 8+< <2 .654 85 E - 8 ) *T *T =^P 5 .’* ! 

V 8»V 8--<X I .5 934 9E-6 ) *P *P/ 5.^ • f Z* J 3 , j 
let V8»VSV<-<S .22ll6E-T3)*T*P*P’»P/ ^ ^ 

t^T V8-V8-<X4-.38031E--^-T.T-T.P/ 59 . 31 ^^^^ 
let V8»V8--<‘a.6l674£-9)»S-S*P*P/C98l3 2)) 

let V 8-V 8 V<-< 9 .6 3403E -ii ) -T -T -» • < I 300) ) 

let V8»V3+'a4.3S64£-6 -5. *P/ C 9 . 31 ) ) 

let V 8-V 3-rC3 .403 97E -4 ) -T -S *P/ C98l 0 ) ) 

f TT V-V4 ♦vs ••■V6+V7 ♦vs 

let V»V»C 130 ) 

OENSITYCS *T^P ) — E0CS>T/?) 

2t R0*CC2.996E^9^P5/ca*996£^9Vl.3l36£^6))t.l368 

LET R0-R0-. 998492- th(P CS-. 697363) 

LET R0»R0-.999973-. 999875 ' - , = . 5 ^^T-T -C6 .79E -3 ) -T-T -T ) 

LEr R0-R0/ CL-C6 .427E-5 )-T;^( 8.0053^-6 )-T-T CO . /yt. o 

INDEXCS^T/P^D— NCS^T>?/L) 
let N-L-L- .32427/ CL*L-C9S2.68t2)‘) 

LET N-N-K.00000L V-C20-.06-T)-S _ 

let N»N-<H0 *972)-' £X? C C -7 .09E-5 ) -T ) 

let n-n+v 


^ derivative 3UBR0UTLME 

NEEDS XU -3) AND Y d -5) 

SLOPE RETURNED AS 19 

let 13* 3 
let 14- a 
LET as- ,0 
LET 16- a 
FOR a*l TO 5 
LET 23*13 ♦XCai 
let Ia- 14+YCII 
LET 15»1S^XC1) •YCi: 

LET 16 -16 ♦XCl 3 -XCa 1 

LET 19* CS -IS -13-14 )/ CS -16-13 -13 ) 

^^JRN EVALUATES D3-CDN/0S) AT CONS i AN i i 

LET YC 3I-N 

let XC 01-S 
LET XC6I-R0 


ORIGINAL PAGE IS 
OF POOR QUALITY 


8915 
8 920 
8922. 
3925 
3930 
8935 
8940 
8945 
8950 
3955 
8960 
8965 
8990 
9000 
9005 
9013 
9012 
9015 
9020 
9022 
9025 
9030 
9035 
9040 
9045 
9050 
9055 
9060 
9065 
9090 


FOR a«l TO 5 

LST S*XC aixa-3 )*<»032} 
GOSU3 8600 
GOSUB 8700 

let xcai*s 
let YCai»N 
NEXT a 
GOSUB 8800 

let D3»a9 
let N»YC 0] 
let 3»XC 0] 
let R0«XC63 

RETURN 
REM 

YC 03«N 
XL 01*P 
XC6 3*R0 
TO 


EVALUATES D2*<DN/DP) AT CONSTANT k.>S 


LET 

LET 

LET 

FOR 


let P=x: 3 3 + <a-3 )»C15 -*‘o 3 
GOSUB 8600 
GOSUB 3700 

let xca3*p 

LET YCai»N 
NEXT a 
GOSUB 8800 
let D2«a.9 
let N»YC 0 1 
let ?»XIL 01 
LET R0*XC61 
RETURN 


9100 

REM 

9105 

let YC 31*N 

91 10 

LET XC 01=T 

91 1 2 

LET XC61*R0 

9115 

FOR TO 5 

9120 

let T»XC 3 1 

9122 

GOSUB 5600 

9125 

GOSUB 3730 

9130 

let xcai»T 

9133 

let YCai»N 

9140 

NEXT a 

9145 

GOSUB. 8800 

9150 

LET Dl«a9 

9155 

let N«YC 01. 

9160 

let T*XC 01 

9163 

LET R0«XC61 

9190 

RETURN 

9200 

REM 

9205 

let YC 3 1»R3 

9210 

LET XC 01=»T 

9213 

FOR 1=1 TO 3 

9220 

let't*xc a 


EVALUATES DL*<DN/DT) AT CONSTANT ?>S 


■►(a-3 )«c 1 ) 


:ES Al *-<ORa/DT )/R0 
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9225 
9230 
9235 
9240 
9245 
9250 
9255 
9260 
9290 
9300 
9305 
9310 
9315 
9320 
9325 
9330 
9335 
9340 
9345 
9350 
9355 
9360 
9390 
9400 
9405 
9410 
94 L5 
9420 
9430 
9440 
9450 
9500 
9505 
9510 
9515 
9520 
95 25 
9530 
9999 


EVALUATES 31 » (UR0/DP )/E0 


GQSU3 8600 
let xcaa»r 
let YLa3-P.0 
NEXT Q. 

GO SUB 8800 

let 

let R0“YC 01 
let T»XC 03 
HETUP-N 
REM 

let YC. 31«R0 
let XC 03»P 

“OR 2*1 TO 5 _ 

let P»Xt 0 3-^<2-3>*5li;^6 ) 
GO SUB 3600 

let xcai«p 
let YCQ3»R0 

NE>T 2 
GOSUB 3600 
let 31*29/ R0 

let R0=YC'0 3 
let ?*XC 31 

return 

REM 
LET 
LET 

let 

LET 
LET 

RETURN 
LET Ul» a 
FOR 23»L TO 3 

let V2»U2-^<33-2)»FI. 
GOSUB 9400 
LET U1 »U1 •►VI 
NEXT 23 ' 

RETURN 
END 


CALCULATES INTENS ITY CV 1 ) 
V1»CVI*U1 )-*-<4*(V2-Hl )*CV2-ht ) 5 
V3*XV1 *Wl )-f'<4*'<V2-*-Hl ) *CV2+Hl ) ) 

V L -VI t^A^VC -A3 *T6 -W8 -Is / C W3-W 8 +4 -V 2 -V2 ) 
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